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o": 

C " 3 , Abstract. The aim of this paper is to examine some obtained exact solutions of 

f^*^ , the Einstein-Maxwell equations, especially their properties from a chronological 

^O ■ point of view. Each our spacetime is stationary cylindrically symmetric and it is 

filled up with an perfect fluid that is electrically charged. There are two classes of 
t-*H' solutions and examples of each of them are investigated. We give examples of the 

QJ ' first class both for the vanishing as well as for the non-vanishing Lorentz force. 

1. Introduction 

>; 

Thanks to its remarkable properties the Godel solution (Q, |g|, Q) became subject 
;^ ' to several interesting generalizations. Banerji and Banerjee in [1 have found solutions 

of the Einstein-Maxwell equations for a charged perfect fluid. Also generalization to 
the case with or without the electromagnetic field and when non-constant component 
9zz or g u of the metric tensor is admited is possible (0 and references citated therein). 
In this paper we give a number of solutions of the Einstein-Maxwell equations with 
the vanishing cosmological constant which are generalizations of the metrics discussed 
in p| and which involve electromagnetic field as well as the non-constancy of g zz or 

$H . gu- These solutions result from two classes of exact solutions of the Einstein-Maxwell 

equations for a spacetime filled up with the charged perfect fluid. The subject of this 

J> , work directly generalizes the paper given by Mitskievic and Tsalakou B, who used 

the H-M conjecture (||), and therefore their designations are partly preserved. 

The plan of the paper is as follows. In section g we obtain the Einstein- 
Maxwell equations for the stationary and cylindrically symmetric spacetime with 



- * — i 



vector potential ( |2.3| ) . Through the whole paper we systematically use an orthonormal 
basis in which the stress-energy tensor of the perfect fluid is diagonal. In section [| 
we get a general solution of the first class that corresponds to the non-constant 
component g zz of the metric tensor. In the case of vanishing Lorentz force we write 
down the inequalities implied by the energy conditions as well as the values of the 
Riemann tensor for this solution. Section H derives a second class corresponding to 
the non-constancy of gu component of the metric tensor. 



Cylindrically symmetric perfect-fluid universes 2 

2. The Einstein-Maxwell equations 

In this paper, we consider a spacetime filled up with the charged perfect fluid of 
the pressure p and the energy density fi. Let us choose local coordinate system 
x^ = (t,ip,z,r) | of the comoving coordinates, where p is angular, r radial and 
z ordinary Cartesian coordinates. We search for the metric which is stationary and 
cylindrically symmetric and depends on the only coordinate r, so we will restrict our 
attention to the case when the spacetime admits Killing vector fields in remaining 
directions, i.e. in t, ip and z directions. Since we demand the metric to be invariant 
over simultaneous time reversion t — ► — t and reflection tp — » — (p, and over the inversion 
z — > — z, the spacetime metric tensor has the form [[]] 

ds 2 = e 2a (dt + f d(p) 2 - I 2 dp 2 - e 27 dz 2 - e 2S dr 2 , (2.1) 

where all metric functions as well as the pressure and the charge density depend only 
on the coordinate r. The function 5 can achieve any value. The velocity vector is 
given by u = e~ a dt and the acceleration of the particles of the fluid is 

u = y u u = a'e~ 2S d 3 , (2.2) 

where comma means derivative with respect to r. The expansion tensor as well the 



shear tensor of the spacetime with the metric (2.1) is vanishing 



S,i*, = ~ £ h M „ = - £ g^ - u (ai u„) = . 

Z u Z u 

Finally, let the one- form of the vector potential be expressed like (see p| ) 

A = m(r) dip + n(r) dt . (2.3) 

Let us introduce an orthonormal basis with the basic one-forms as follows 

e 6 = c a (dt + /d^) , e i = id(p, 

6 2 = e 7 dz , 9 s = e s dr . 



(2.4) 



Using (2.2) for mixed components of the two-form of the electromagnetic field 
one gets 

F; = l~% { [n'l 2 e- 2a + (m' - fn')f] % (m' /n')^}+e- 2 "^(«'<5°+m'^).(2.5) 

The Maxwell equations with sources 5F — 4irj give 

l- 1 e- a -' f - s {e a+ ^ s r 1 [[n'l 2 e~ 2a + (to' - fn')f\5% - (m 1 - fn')d%]}, 3 = -Airf, (2.6) 

j(r) being the current density, Sj = 0. We postulate that the fluid particles carry the 
charge, j — pu, p is the charge density. It will be accomplished if 

(m' - fn')c a+ '<- & = Bl , (2.7) 



B is constant. The equations ( |2.6| ) and (2/7) result in the following expression for the 
charge density 

m — jn dr \m — fn J 



The choice (2.3) corresponds to the formulae for the electric and the magnetic 
fields (|) 

E f = n'c- a - s , B 5 = 5e- Q ~ 7 . (2.9) 

§ The coordinates will be numbered (t, ip, z, r) = (0, 1, 2, 3). 

|| We use the units with the speed of light c and the Newtonian gravity constant G equaled to one. 
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With the help of the constraint (^J) the vorticity one-form is given by 

■dz . 



1 / , x Bf 

2* {uAdu)= 2W~fn<) 



The total stress-energy tensor of the spacetime T to tai arises as a sum of the stress- 
energy tensor of the electromagnetic field 

Teimag = (8tt) ~ l l~ 2 &- 2S { 2n' (m' - fn')le~ a {Q b ® 9 1 + 9 1 ® 9 5 ) 

+ [n' 2 Z 2 e" 2a + (m' - fn') 2 ] 9 fi ® 9 fi + [n' 2 ; 2 C - 2 « + (m' - /n') 2 ] 9 1 

+ [n ,2 l 2 e- 2a - (m' - fn'f] 9 2 ® 9 2 - [n' 2 ; 2 e~ 2a - (m' - /n') 2 ] 9 3 
and the stress-energy tensor of the perfect fluid 

Tfiuid = /i9 6 <g> 9° + p (V <g> 9 1 + 9 2 <g> 9 2 + 9 s <g> 9 s ) . 
The Bianchi identities reduce to the single equation of the motion of the fluid 

p' + a'(p + fx)+pn' = 0. (2.12) 



(2.10) 
9 1 
9 s } 

(2.11) 



The non-zero components of the Einstein tensor in the orthonormal basis (2.4) 
can be written in the form 



Wc£ — "i0 



_ \ -2a-y-S _^_ (r-lfl a 2u ,,<<--! ^. 

2 dr V ■* 



c 



G^ - G ff = ie—T-'^- (r^Q + 7 )'e Q+ ^) - 2a' 7 'e- 25 



G« - 



" G;g,g 

G" : 



-ig-" — y-<5 



; _d_ 

dr 



dr 



(' 2 



2 dr 



r l e a- 7 -i - ( Z 2 e 2-y) 



dr 



+ I r 2 / /2 e 2(a ! -5) 



G«. = (a'/' + i'V + h'a') rh- 25 + Ir 2 /' 2 e 2 ( Q -^ , 
and the Einstein-Maxwell equations G^c — QtiT^o are 



_d_ 

dr 

_d_ 

dr 

l__d_ 

2 dr 

l_d_ 

2"d7 



Jla fl 



r 



ml — frt! 
(a + 7 )' 



fn> 

-la 



- fn' dr 
^ d_ 

- /n' dr 



An' , 

2a / 7 / 
m' — fn' 

' (Z 2 e 2a ) 
(Z 2 e 2 ^) 



2n' 2 e 



/2„-2a 



?n' 



2A" 



/n' 

167d 2 pe 



m! — fn 1 

8Tr(p - M) j2 e 2* 



1 /' 2 e 2 



/n' 



2 m' — /n' 



2n' 2 ; 2 e" 2Q 



/n' 



m 



(a' + 7')^r- + 2a' 7 '/ 2 = 167r/ 2 pe 2<5 
dr 



-/' 2 e 2a - 2n' 2 l 2 c- 2a + 2(m' - /n') 2 



(2.13a) 
(2.136) 
(2.13c) 
(2.13d) 
(2.13e) 



In summary, we obtained the system of six equations (2.13a)-(2.13e) a nd ( 2.6) 
minus (|2.13e ), which is an integral of the (2.12). After eliminating of p from ( 2.13c| )- 
( [2.13c} ) with the help of the (2.13e) one gets the system of six equations for totally nine 
unknown functions a, 7, 5, f, I, /1, p, m and n. It means that three out of nine functions 
can be chosen arbitrarily. These three degrees of freedom are equivalent to introducing 
S, p and nip). In the next we will consider the m and n fixed, and moreover we impose, 
as in H, aditional condition: either a or 7 are constant, which is in accordance with 
two classes of solutions. 
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3. First class a is constant 

When a — const the Einstein-Maxwell equations ( 2.13c) )-( 2.13cj ), ( |2.6| ) can be easily 
integrated in terms of the electric and the magnetic potential to and n. Equation (2.12) 
now says that the Lorentz force is balanced by the pressure gradient, gradp + pE = 0, 
resulting in the geodesic motion of the fluid particles. We can, without loss of 
generality, suppose that a is zero. The metric components are 



7 



2n- A 



to ne 
ii 



2,i 



dr + F 



2n' 



dr + C 



(m! — fn') dr , 



m' — fn' 
l 2 e- 2 ~< = E-A f (m-fn + k)(m'~ fn') e - 2 ' f dr. 



(3. 


la) 


(3. 


16) 


(3. 


lc) 



I n (3.1c ) C, E, F, k are constants of integration. Inserting ( 3.1a ) into ( 2.13e ) and 
( ^.13 4 ) yields for the pressure and the energy density 



8irp = B 2 



8?r/x = B 2 



An 2 



Y 2 + n' 2 p _ 2 , m- fn + k _ ^ 



An 2 - 37 



(to' — fn') 
a , 5n ,2 



(to' — fn') 



rl 2 + 6 7 



to' — fn' 

,to — fn + k 
to' — fn' 



a -27 



-2 7 



(3.2) 



3.1. Example of the first class with non-vanishing Lorentz force 

We give here explicitly the solution when to = rn + %[3n 3 , r and (3 being constants, 
and k, F and C in ( 3.1c| )-(3.1c) are zero. In this case (according to the (2.5) the 
Lorentz force does not vanish) the metric can be written like this 



ds 2 



~dt+ (2f3n 2 



0) dip] 1 



il 2 



(3£/rV" 2 



An 2 



1 d^ 2 



e 2 " dz 2 



3n' 2 e 2n dr 2 



B 2 3Ef3- 2 e 2n2 -An 2 + 1 ' 



(3.3) 



and for the charge density (2.8), for the the pressure and the energy density (3.2) one 
has 



8irp = B 2 
8nn = B 2 



E 



An 2 



i)- r 



-2n- 



np = —B n 



1 



-2n" 



{[f" 2 -i]- 2 -*-|( i2 » 2+5 >} 



(3.4) 



3. 2. Solution with the vanishing Lorentz force 



Metric given by (3.1o)-(3.1c) is not too transparent. In the rest of this section we 
restrict ourselves to purely magnetic field when the electric potential n is constant, 
which will be marked as jb, i.e. when according to ( |2.5| ) or (2.7) the Lorenz 
force F^u u vanishes. The result written in terms of the magnetic potential to is 
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^(F was gauged away) 

ds 2 = ( dt + bm dip) 2 - {Ec 2Cm + Am + v) dip 2 - c 2Cm dz 2 



1 me 



l2 2Cm 



dr 2 , (3.5) 



(3.6) 



B 2 Ee 2Cm + Am + z/ 
A-b 2 _ A + Ak 

The charge density, pressure and the energy density are 

Airp = -B 2 be- 2Cm , 8irp = B 2 C 2 E , 

8irfj, = B 2 [{b 2 - 2)e- 2Cm - 3C 2 E] . 

Our fluid moves geodesically, without shear and expansion [j. For m being 
increasing or decreasing function defined on (0, oo) and covering the whole range 
(0, oo), it can be treated as a new radial variable and one obtains the first family of the 
solutions in 14] , if we interpret the coordinate ip as an ordinary Cartesian coordinate 
(i.e. one abandons the periodicity of the ip coordinate). The same will be true in the 
section fl The solution without the electromagnetic field (B = 0) was recovered by 
Wright (1). 

3.2.1. Energy conditions. Dominant energy condition implies inequalities 

{b 2 - 2)c~ 2Cm - 2C 2 E > , (3.7a) 

(6 2 - 2)c~ 2Cm - AC 2 E > . (3.76) 

In our case the strong energy condition is satisfied if the dominant energy condition 
is. Altogether the energy conditions will be fulfilcd if (3.7c|) or (3^73) is satisfied 



depending on whether E is negative or positive respectively, or b > 2, if E is zero. 

3.2.2. Curvature. For the future convenience we mention the concrete values of the 
Riemann curvature tensor for the metric ( |3.5| ). All non- vanishing components of the 
Riemann tensor are given by 

Rfzfz = R$m = \B 2 e~ 2a (XCe- 2Cm + 2C 2 E) , 



.,., r .---tf~/r<- v 



r> _ p _ R 2.2„-4a„-2Cm 

Htt.u. — nt-„;~„ — — — n o c 



Rir$r = -Ri^ z = lB 2 bCe- 3a (Ee 2Cm + Am + ^e~ 2Cm 
R^ f = - -B 2 c- 2a (3b 2 e- 2a + 2\C)e~ 2Cm + B 2 C 2 Ee~ 2a 



(3.8) 



If In fact, there is another way of deriving localy the same metric as (3.5) that introduces new 
variable by rescalling the coordinate r by the definition m(r) = m. Then for the two form of the 
electromagnetic field F = dm A dip. In effect, this substitution r esul ts in inserting unity except m' . 
+ Thanks to the homogeneity of the pressure ( p.q ),the solution ( p.5| ) can be reinterpreted also as a 
solution describing charged dust in the spacetime with the cosmological constant A = —B 2 C 2 E (the 
Einstein equations in our convention are G — Ag = 8ttT). 
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3.2.3. Closed timelike curves. For the analysis of potential chronology violation we 
use the Carter theorem (Q, ||). In our case (3.5) the isometry group, generated by the 
Killing vector fields d t , d v , d z , is Abelian and timewise orthogonally transitive, since 
the hypersurfaces of transitivity given by constant r are everywhere timelike, g rr < 0. 
The only exception could be the case C = 0, but then the chronological structure 



follows from the continuity and we will meet such an example in subsubsection 3.2.4 



Chronology will be preserved when one is able to find such constants p, q, s that the 
differential form ijj = pdt + qdip + sdz is everywhere timelike or null. Since the 
contribution of s is always negative, it can be treated as zero. 



3.2.4- Example of the first class with vanishing Lor entz force. Let the vector pote ntia l 
be A — 2a 2 Bsh 2 (r/2a) dip. Here a is a length characteristic^], B is according to (2.9) 
a value of the magnetic field on the rotation axis. The choice of the constant of 



integration E = -v = (2a 2 B 2 C 2 )' 1 
2a (Aa 2 B 2 + 2 - 2C) k 



tf = 4 - 4CB- 



ds 2 



dt 



e 2Cm dz 2 



2a 2 C 2 m-e 



mdip 

'2,;2Cm AJ2. 



2a 2 r 



2Cm 



2(aB)~ 
-2C(1 



yields the metric 



C)m - 1] dip 2 



o2Cm 



1 



(3.9) 



2C(1 - C)m 

where, for convenience, we have denoted 2a 2 BC — » C and m = sh (r/2o). For the 
physical quantities fi, p and p one obtains 

1 _ @\ p -2Cm ^_ 

2a 2 



87r/i = 2[B 



2irp 



-B 



B 2 



\-C 



(3.10) 



-2Cm 



8np = 



2a 2 J ' r 2a? ' 

Signature in (|3.9| ) is always correct and the energy condition ( |3.7&|) requires C < 0. 
From equations (3J5) we can see that for negative C the physical singularity occurs 
when r — > oo. From (|3.9| ) one can also see that no event horizon is present for any 
finite value of r. In the limit C — we get Bancrji-Banerjee solution (formulae (11) 
and (12) in the alternative interpretation) in |3J which is spacetime homogeneous and 
singularity free. So we can interpret C as an indicator of the difference of the spacetime 
from the spacetime homogeneity. 

In fact, the choice m = sh {r /2a) represents the family of models in which m is 
increasing or decreasing function with domain (0,oo) and image (0, oo). Each such 
choice after the introducting of the new variable m = m(r) leads to the (3.9). From 
this point of view it is obvious that for example the choice m oc r -1 has the same 
physical content as m = sh (r/2a). The same remark will be true in subsection 4.3. 

The condition V'qV'" > of the subsubsection |3.2.3 reads 

-\ 2 



1 



— ie 2 ^~2C(l-C)m-l}- 



(Aa 2 B 2 + 2 - 2C) 2 m - — 



2ap 



>0 . 



(3.11) 



Because q must be zero in order not to predominate in (3.11) for small r, the condition 
for the non-existence of the closed timelike curves (hereafter CTC) has the familiar 
form g w < (or g lpv < for the non-existence of the closed causal curves), and we 
see that CTC will apper always for sufficiently large r. 



With respect to the previous footnote, 



(2A) holds in the alternative reinterpretation. 
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Table 1. Numerical dates that relate the quantity of the magnetic field B on 
the rotation axis to the radius R of the first null circle (£, z, r = const), and 
to the cosmological constant A in the alternative interpretation. For simplicity, 
dates are given only for C = 0. The matter density used is 10 -26 kg ■ m — 3 , which 



corresponds to the period of the rotation 7 • 10 10 years. 




B [10 3 Gauss] 0,0 3,7 4,8 6,7 8,2 


10,5 


A [10- 53 m" 2 ] -9,31 -8,20 -7,45 -5,59 -3,73 
R [10 8 ly] 137 132 129 123 118 


-0,28 
110 



Thanks to absence of any event horizons, for any two spacetime points p and q , 
one has p » q and simultanously q S> p (@)- Especially I + (p) = I~(p) = J + {p) = 
J~{p) — M. Although through every point passes some CTC, each such curve must 



inevitable cross the region r > R, where R is determined as a root of (3.11). The CTC 
we are interested in are nontrivial (|§). Moreover, because the Carter theorem can be 
applied independently on whether or not tp is periodic, one can see that CTC occur 
even if we interpret tp as the ordinary Cartesian coordinate. (The chronology could be 



preserved in the case that we took only part of the (3.9) confined to the region r < R 
and tried to match it to some other chronologically well behaved solution.) 

The magnetic field has always positive influence on the chronology violation, 
precisely the larger value of the magnetic field on the rotation axis the closer to it will 
be the "first" null circle followed by CTC (see table table |l|). 

Similarly, another Banerjee-Banerji solution (expression (15) in p|) or the Som 
and Raychaudhuri solution O] can be got by suitable choice of the constant of 
integration. 

4. Second class of the solutions 7 is constant 

Solution with constant 7 will be refered as second class and here we restricted ourselves 
only to the case when the Lorentz force vanishes, with the electric potential n — 6/4. 



Equation (3.1c) remains valid in this case too. 7 can be put to zero and the solution of 
the Einstein-Maxwell equations (2T3q)-(2T3q) and (2.6) can be written in the form 



/ = -— c- 2Cm + H ,c 2a = e 2Cm , 

I? (4-1) 

l 2 = ^e- 2Cm -2m 2 + Dm + E, 

with constants of integration C, H and D = Hb — 4k. For the energy density, the 
pressure and the charge density one has 

16tt^ = B 2 e~ 2Cm [CD + 2(1 - 2Cm)] , 

16irp = B 2 e- 2Cm [CD - 2(1 + 2Cm)] , Anp = -B 2 be~ 3Cm . 

Thanks to inhomogeneity of the pressure (§^3), the fluid already does not move along 



geodesic lines. Its acceleration (2.2) is 

V u u = Cm'e- 25 d 3 . 

In contrary with the corresponding solution of the first class ( |3.9| ) , because of the 
inhomogeneity of the pressure, the metrics of the second class cannot be reinterpreted 
as describing dust with the non-zero cosmological constant. 



Cylindrically symmetric perfect-fluid universes 
4-1- Energy conditions 



In the case of the second class (4.1) and (4.2) the energy conditions altogether will be 
fulfilled if following inequality will hold 

CD > ACm . (4.3) 

4-2. Curvature 

The non-zero components of the Riemann tensor are given by 



R ifif = R i<pi0 - b2Cc 



2m e 

2 



2Cm 






C [ 2m - ~ 



(4.4) 



-2Cm 



4-3. Example of the second class of the solutions 

This special solution corresponds to that discuss ed in |3.2.4 in the sense that when C 



goes to zero, we obtain the same solution as in |3.2.4| when C goes to zero, i.e. (11) 
and (12) in ||. For let us choose A = 2a 2 i?sh (r/2a) dip and constants of integration 
b 2 = 4 + 2 {aBy 2 , D = 2B- 1 + 2C*" 1 + (a 2 B 2 C) _1 ,E = -b 2 (2C)~ 2 ,H = b (2Cy 1 . 
Components of the metric tensor read (denoting 2a 2 BC — ► C and m = sh 2 (r/2a)) 



I 2 = 4a 2 



2a 2 B 2 



,2<< 



1 



,2Cm 



, / 



(4a 2 B 2 



2Cm\ 



- 2a z B 2 m 2 



C 
(2a 2 B 2 



(1-e 

f 1 + C)m 



2a 2 B 2 



1 



(4.5) 



2C 2 ' C 2C 2 

The energy density, the pressure and the charge density are given respectively by 



8^ = ( 2B A + 

l + C 



1 + C 

2a 2 



2B 2 Cm ) e 



-2Cm 



8irp = 
2np = 



2a 2 



- 2B 2 Cr 



-BIB 2 



1 
2^2 



-2Cm 



-3Cm 



(4.6) 



It can be seen from the form of metrics fl4.5| ) that the signature is correct for 
C < and, in case when B is zero, then for every C. The absence of the event 



horizons for any finite value of r is apparent from the form of the metric (4.5). From 
analysis of the energy conditions ( |4.3| ) it follows that these will be fulfilled for intervals 
— (2a?B 2 + 1) < C < 0, and if the electromagnetic field vanishes for C > — 1. One 
can convince oneself that CTC will occur if — (2a 2 B 2 + 1) < C < when B =/= 0, 
and C > — 1 when B = 0. Every two points of the spacetime can be connected each 
other by both future and past directed timelike curve which is non-trivial, so that the 
spacetime is totally vicious (|§). Formulae (4.4) show that the physical singularity 



occurs for C < when r tends to the infinity. The only exception is the case with the 
vanishing electromagnetic field (B = 0) and C = — 1, which is after transformation 
e -2m _ 1 — ^2^2^ w j^ n (jj- 2 = 2a 2 , the Minkowski spacetime in the rotating cylindrical 
coordinates ip — > ip — cut. 
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5. Conclusion 

From the requirement of the stationarity and cylindrical symmetry and from the 
assumption that the spacetime is filled up with an perfect fluid that is charged, we 
have obtained the system of ordinary differential equations (2.13tj)-(p.f3c() plus (p.6|). 



that was solved in terms of the functions m and n in two cases: the first class (|3.f (\) 



([j.f (j), (|3.2[) with constant a and the second class (4.f ), (4.2) with constant 7 (but 



we gave the solution of the second class only for the vanishing Lorentz force) . These 
two general classes are both shear-free. Fluid of the first class moves geodesically 
and its velocity is the Killing vector field, ft is not true in the second class thanks 



to the inhomogeneity of the pressure (4.2). We shown explicit special solution of the 



first class for the non- vanishing Lorentz force (formulae ( p.3[ ) and (3T)), but without 



detailed analysis, and in some details we discussed the example of the first class 



(su bsubsection 3.2.4 ) and the corresponding example of the second class (subsection 
4.3) with the vanishing Lorentz force, that are generalizations to the solution given 
in jj|. These metrics with three parameters (the length characteristic a, the value of 
the magnetic field on the rotation axis B, and C determining the difference of the 
solutions from the homogeneity) contain CTC at least for some interval of the values 
of C ' . ft turns out that the magnetic field has positive influence on the appearance of 
the non-trivial chronology violation (in fact by sufficiently large magnetic field we can 
always ensure the chronology violating) . 
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